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Multi-strange baryon production in pp collisions at y/s = 7 TeV with 



A measurement of the multi-strange E~ and £1 baryons and their antiparticles by the ALICE experi- 
ment at the CERN Large Hadron Collider (LHC) is presented for inelastic proton-proton collisions at 
a centre-of-mass energy of 7 TeV. The transverse momentum (pj) distributions were studied at mid- 
rapidity (|y| < 0.5) in the range of 0.6 < pj < 8.5 GeV/c for E~ and baryons, and in the range of 
0.8 < pj < 5 GeV/c for Q. and Q. . Baryons and antibaryons were measured as separate particles 
and we find that the baryon to antibaryon ratio of both particle species is consistent with unity over 
the entire range of the measurement. The statistical precision of the current data has allowed us to 
measure a difference between the mean pj of E~ and £2 (£2 + ). Particle yields, mean pj, and 
the spectra in the intermediate pj range are not well described by the PYTHIA Perugia 2011 tune 
Monte Carlo event generator, which has been tuned to reproduce the early LHC data. The discrep- 
ancy is largest for £2 This PYTHIA tune approaches the pj spectra of E~ and E + baryons 
below pj < 0.85 GeV/c and describes the E~ and E + spectra above pj > 6.0 GeV/c. We also illus- 
trate the difference between the experimental data and model by comparing the corresponding ratios 
of (i2~+£2 + )/ as a function of transverse mass. 
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1 Introduction 

The multi-strange baryons, Q. (sss) and S~ (dss), are particularly important in high energy particle 
and nuclear physics due to their dominant strange quark (s-quark) content. The initial state colliding 
projectiles contain no strange valence quark, therefore all particles with non-zero strangeness quantum 
number are created in the course of the collision. Moreover, the energy of the Large Hadron Collider 
(LHC) and its high luminosity allow for an abundant production of strange hadrons. These two factors 
make multi-strange baryons a valuable probe in understanding particle production mechanisms in high 
energy collisions. 

We present a measurement of Q. and Q. + baryon transverse momentum (pj) spectra and yields in proton- 
proton (pp) collisions at a centre-of-mass energy (y/s) of 7 TeV, a measurement of S~ and S + yields and 
spectra at the same energy, and a comparison of these data to a recent pp event generator, PYTHIA 
Perugia 2011 central tune (P2011). The measurements were obtained using the ALICE experiment fl] 
at the LHC. 

ALICE is a general purpose detector designed to study both pp and Pb-Pb collisions at TeV-scale ener- 
gies. A six-layer silicon inner tracking system (ITS) and a large-volume time projection chamber (TPC) 
enable charged particle reconstruction with excellent momentum and spatial resolution in full azimuth 
down to pj of 100 MeV/c. 

2 Data sample and cascade reconstruction 

Multi-strange baryons are studied in a sample of approximately 130 million minimum bias y/s = 7 TeV 
pp events, collected during the 2010 data taking. The sample is corrected for trigger inefficiencies and 
biases to recover a normalized sample of inelastic (INEL) events, as described in The events are 
selected within 10 cm of the detector's centre along the beam direction, with vertex resolution in the 
transverse plane of a few hundred micrometres. The event vertex range is selected to maximize particle 
trajectory (track) reconstruction efficiency within the ITS and TPC volume. 

S~ and S + (S 1 * 1 ), as well as Q. and Q. + (H ± ) candidates are reconstructed at mid-rapidity (\y\ < 0.5) via 
their characteristic weak decay topology, £~(£ + ) — > A(A) + n~(n + ), and £1 (£1 ) — > A(A) + K~(K + ), 
as described in detail in 0. The branching ratios for these decay channels are 67.8 % for baryons 
and 99.9 % for E ± . Charged particles, compatible with kaon, pion and proton hypotheses, are identified 
using their energy loss in the TPC. The topology of the £2 and S~ weak decay is cascade-like and 
consists of a V-shaped decay of the daughter A baryon (A baryon hypothesis is identified as a "V0") 
plus a negatively charged track (h~). The same applies to antibaryons, however in that case the decay 
products are the A daughter particles and a positively charged track (h + ). In general, the acceptance and 
efficiency depend on both y and pj. We chose the y interval such that our efficiency and acceptance 
depend only on pj. Candidates are selected by placing restrictions on the topology of the decay. These 
have been optimized to obtain maximum mass signal significance and are listed in Table Q] 

The resulting invariant mass distributions for both species hypotheses are shown in Fig. Q] The signal 
extraction method is described in detail in J3]. The signal is extracted using a bin-counting method 
and then corrected for detector efficiency and acceptance using PYTHIA Perugia [4] generated Monte 
Carlo events propagated through ALICE using GEANT3 10. 

3 Systematic uncertainties 

There are two types of systematic uncertainties in the resulting particle spectra: /^-dependent systematic 
uncertainties that are due to the efficiency determination and the signal quality at a given pj, and the pj- 
independent uncertainties due to normalization and other factors explained below. 
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Table 1: Selection criteria parameters for VO (A) and cascades (S and Q. ) presented in this Letter. If a criterion 
for and finding differs, the criterion for il* hypothesis is in parentheses. DCA stands for "distance of 
closest approach," and PV for "primary event vertex." The fiducial volume is defined by the coordinate of the decay 
vertex position, the transverse radius, R2d- 9 is the angle between the momentum vector of the reconstructed VO 
or cascade, and the line segment bound by the decay and primary vertices. For cascades, the curvature of the 
particle's trajectory is neglected. 



VO finding criteria 



DCA (h* to PV) 
DCA (h~ toh + ) 
A mass (myo) 
Fiducial volume (7?2Z>) 
VO pointing angle 



> 0.04 (0.03) cm 
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Fig. 1: The invariant mass distributions of Z~ (a) and Q. (b) baryon candidates (solid histograms) and their 
antiparticles (dashed histograms). Also marked (in shaded blocks) background sampling regions used in the signal 
extraction. The entire measured pj range is presented. 

The point-to-point systematic uncertainties vary between 1-4 % for H ± , and 1-9 % for H ± , with minimum 
uncertainty found at pj = 1.5-4.0 GeV/c for both species. The /?x-independent systematic uncertainties 
stem from several sources and reflect the following: 



- the uncertainty in determination of the material thickness traversed by the particles (material bud- 
get), 4 %; 
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- the use of FLUKA OQ to correct |0 the antiproton absorption cross section in GEANT3 Q, 



- the uncertainty in TPC particle identification via energy loss, 1.5 %; 

- the uncertainty on the track selection in the TPC, through the restriction on the number of TPC 
pad plane clusters used in particle reconstruction, 3 %; 

- in the case of O , the removal of cascades that fit the S baryon hypothesis, 1 %. 

The limited /^-coverage and determination of the total number of inelastic events used for yield nor- 
malization lead to an additional uncertainty in the particle yields and mean pj ((pt)) values. The INEL 
normalization [ 2 ] leads to a +7.0 % and -3.5 % uncertainty on the yield for all measured particles, while 
the limited pj coverage causes a 4.5 % uncertainty on the (pj) of all species, 5.5 % uncertainty on 
the yield of E ± baryons, and 6.5 % on the yield of O . While the systematic uncertainties (both pj- 
dependent and /^-independent) associated with each spectrum point affect the determination of (pj), 
the systematic uncertainty on the (pj) for all species is dominated by the 4.5 % error due to the limited 
px coverage. Similarly, the systematic uncertainty on the yields is dominated by the uncertainties due to 
low-px extrapolation and event normalization. 

4 Results 

4.1 Corrected pj spectra and Tsallis fits 

The corrected multi-strange baryon yields per pj bin per unit rapidity (1/Nhmel x d 2 N /dydpj) are shown 
in Fig. 0a). They span from pj = 0.6 to pj = 8.5 GeV/c in the case of S~ and S baryons and from 
pj = 0.8 to pj = 5 GeV/c for Q. and Q + baryons. The Tsallis function is used for fitting the spectra, 
as the measured pj range covers both soft-physics and fragmentation particle production regions. The 
functional form is shown below: 



where T, n, and dNIdy {dN/dy representing the particle yield per unit rapidity) are fit parameters, 



The function is grounded in Tsallis statistics ; it approximates an exponential component (represented 
by the T parameter), as well as a power-law dependence for the high-px tail. In Tabled we list the fit 
results for each particle and antiparticle and the corresponding extrapolated dA7d;y and (pj). 

The central values of the fit parameters, listed in Table |2] are obtained using the statistical error only. The 
\ow-pj extrapolation of the yield from the Tsallis fit is ~ 23 % for E ± and ~ 26 % for Q^. The value 
of (pj) for each particle was computed using the fit over the entire pj range including the extrapolation. 
The antiparticle to particle ratios were found to be compatible with unity at all pj. 

4.2 Excitation functions 

Our measurements of multi-strange baryons can be placed within the broader context of existing pp 
collision data. We compare to multi-strange baryon yields in pp collisions measured by the STAR Col- 
laboration at y/s = 0.2 TeV 1 10]. and also to the data obtained by ALICE and CMS at y/s = 0.9 TeV 
Oil iTTTTl . There are also data from pp collisions, obtained by the CDF lfT2l and UA5 lfl3l Collaborations. 
We omit the comparison to these data due to a significantly different kinematic range of the experiments. 



1 %; 




mj = \ + p\, and mo denotes the particle mass. 
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Fig. 2: (a) E~ and £1 baryon (solid circles and squares, respectively) and their antiparticle (open symbols) spectra, 
shown with Tsallis fits, (b) Experimental data to Monte Carlo (PYTHIA Perugia 201 1) comparison. The errors are 
added in quadrature. The normalization uncertainty is shown as a black band. 



Table 2: Tsallis fit parameters, yields, and (pj) for each particle species at |y| < 0.5, as well as yields and 
(pj) extracted from PYTHIA Perugia 2011 2] simulations. Statistical and then systematic uncertainties for the 
experimental values are listed. Sufficient statistics were generated to have less than 1 % error on all model values. 



Particle 


T (MeV) 


n 


X 2 /NDF 


ANI&y x 10 3 
data 


<jn) (GeV/c) 
data 


dN/dy x 10 3 
P2011 


(p T ) (GeV/c) 
P2011 




344±5± 10 


10.8±0.4±0.8 


17.4/15 


8.0 ±0.lj#5 


1.21 ±0.01 ±0.06 


5.38 


1.02 




339±5±9 


10.4±0.4±0.5 


14.4/15 


7.8 ±0.1+^ 


1.21±0.01±0.06 


5.21 


1.02 




460 ±40 ±60 


20±9±8 


8.8/5 


0.67 ±0.03+°$ 


1.47 ±0.03 ±0.09 


0.276 


1.14 




430 ±30 ±40 


14±5±6 


7.0/5 


0.68 ±0.03+°$ 


1.44± 0.03 ±0.08 


0.266 


1.16 



For STAR, ALICE, and CMS data, an increase in dN/dy as a function of collision energy is observed, 
presented in Fig. Ha). We note that the CMS Collaboration used non-single-diffractive events (NSD) to 
normalize the yield, while in ALICE a normalization to the inelastic events (INEL) was used. For a direct 
comparison at LHC energies, the INEL E 1 * 1 yield has to be scaled up by 26 % to get the yield normalized 
to NSD events [2]. After scaling, the S yields per unit of rapidity obtained by ALICE agree with those 
published by CMS ifTTl . For S~ baryons and antibaryons, we also observe a slight rise in mean pi with 
collision energy, as seen in Fig. Hb). The (pj) of fl ± baryons at -y/s = 7 TeV is consistent with 0.2 
TeV data, where O and (pj) were consistent within large experimental error. Due to the precision 
of the current measurements, a significant separation between the (px) of and is observed in 
\fs = 7 TeV pp collisions. 

4.3 (Q. -+Tt) I (E +E + ) ratio 

The composition of S~ and Q. baryons differs only by one valence quark flavour: the <i-quark in S~ is 
replaced by the s-quark in Q. . To investigate possible differences in the production mechanism of multi- 
strange baryons with and without the non-strange quark, we study the ratio of (Q.~+Q. + ) to (S~+E + ) 
baryons as a function of pj. The dependence on particle mass is reduced by constructing spectra as a 
function of (mj — nio) for each baryon species. To increase the statistical significance of the measurement, 
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Fig. 3: (a) dJV/dj and (b) (p T ) of S ± and Or 1 as a function of collision energy. The STAR and CMS data are 
normalized to NSD (see text) events, STAR E^ and il* are represented by open rhombuses and stars, respectively. 
CMS E measurements are shown as open triangles, and ALICE E and as filled circles and squares. Multi- 
strange baryons produced using PYTHIA Perugia 2011 simulation (E* baryons as a long-dashed curve and ii* 
baryons as a dashed curve) are plotted for reference. The uncertainties are added in quadrature. 



£l ± /Z ± 

ratio of fit functions 



Qr/E* PYTHIA P2011 



I— 

E 



0.12 


i i i i i 
- pp,(s = 7 TeV 

7 M<0.5 




0.1 


M 






0.08 
0.06 


y- 




i _ 


0.04 











0.5 1 1.5 



2.5 3 3.5 4 

m T -m (GeV) 



Fig. 4: (Qr+Q. ) to (E~+E ) ratio in *Js = 7 TeV pp events as a function of (mi — mo). Experimental data 
(closed symbols: data points; dashed curve: ratio of Tsallis fits), and PYTHIA Perugia 2011 simulation (solid 
curve). Errors on experimental points were added in quadrature. 



for this ratio, the particle and antiparticle spectra are combined. The ratio of the combined spectra, 
(Q.~+Q. + ) I (S~+E + ), is shown in Fig. |4] We observe an increase in the ratio up to (nix — mo) ~ 1.5 GeV 
(which corresponds roughly to pj of 3 GeV/c for either of the baryons), with a possible slope change at a 
higher (mj — mo). The ratio was composed to investigate the possible saturation of the s-quark production 
with respect to production of the non-strange quarks, which would be indicated by the flattening of 
the (£2~+£2 ) / (S~+S + ) ratio. We found that the currently available data does not allow for a firm 
conclusion. 
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4.4 Comparisons to PYTHIA Perugia 2011 

The production of strangeness in pp collisions is not well described by the currently available models. In 
particular, we compare the obtained data to particle spectra from PYTHIA lfl4l . an event generator based 
on the leading order (LO) perturbative Quantum Chromo-Dynamics (pQCD). PYTHIA is available in 
different tunes, for example those listed in ||4], each reflecting a distinct aspect of particle production 
inferred from experimental data. Several tunes were tested, among them PYTHIA Zl, Z2 lfT31l . and 
Perugia [4] tunes. These tunes were several times to an order of magnitude below the measured multi- 
strange spectra and yields (up to a factor 4 for H ± , 15 for Q^) and thus were abandoned. 

In this Letter, we use the central PYTHIA Perugia 2011 (P2011), one of the more recent PYTHIA 6.4 
tunes, which describes the 7 TeV pp charged particle spectra reasonably well. P201 1 is tuned to the mul- 
tiplicity and charged particle pj distributions from 2010 LHC data (as were other tunes tested), utilizes 
the CTEQ5L partem distribution function, and differs from other PYTHIA tunes by a significant increase 
in multi-strange baryon yields. This is achieved mainly by removing the baryon suppression inherently 
present in the built-in "pop corn" meson creation mechanism, but also by tuning the relative u and d 
vs. s-quark production rates, and adjusting the suppression of the diquark-antiquark hadron production 
scale [4]. The "pop corn" mechanism, used to describe the hadron production in e + e~ collisions via 
chromoelectric flux tubes lTT6l . suppresses baryon production by favouring soft quark-antiquark pairing 
into mesons ifTTl . The mechanism was removed in order to approximate more closely the multi-strange 
yields from LEP experiments |@]. 

Although the charged-particle multiplicities are reasonably described, as mentioned above, PYTHIA 
tunes tend to be several times to an order of magnitude below measured multi-strange values. P2011 
significantly underestimates multi-strange particle yields, as seen in Table |2j and does not reproduce the 
spectral shapes of either S~ or Q. baryons, with two exceptions. The model describes the high pj tail 
of the S 1 * 1 distribution and approaches the S 1 * 1 distribution below pj < 0.85 GeV/c, as shown in Fig. Efb). 

P201 1 also underpredicts (pj) of multi-strange baryons at all energies (Fig. |3jb)), and incorrectly models 
the increase in dN/dy as a function of centre-of-mass energy (Fig. Ha)). Indeed, in experimental data 
dN/dy increases by nearly a factor of three from yfs = 0.2 TeV collisions to those at yfs = 7 TeV 
(factor 35 increase in energy), while P201 1 predicts a more modest gain. In both experimental data and 
P2011, a power-law increase of S~ baryon yield is seen as a function of y/s. Moreover, P2011 does 
not reproduce the relative Il ± /E ± spectral shape, nor the absolute value, although the ratio does increase 
with increased pj, as shown in Fig. |4] 

5 Conclusions 

Our precise measurements of S - , E + , D. , and £2in y 7 ! = 7 TeV pp collisions are a benchmark for 
improving future modelling efforts, including valuable checks on possible hadron production mecha- 
nisms, such as the flux-tube mechanism. In addition, the pj reach of the data to model comparison is the 
highest ever achieved for multi-strange baryons. The relative production of doubly-strange vs. triply- 
strange baryons introduces a further constraint on the pj dependence of particle production from flavour- 
differentiated quarks. These considerations may enable a better insight into pp collision dynamics, which 
in turn will serve as a reference for better understanding of fundamental interactions underlying particle 
creation mechanisms in pp collisions. 
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